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ABSTRACT 


A  three-level  quasi-geostrophic  atmospheric  model  on  a 
mid-latitude  3-plane  with  a  single  wave  in  the  x-direction 
and  a  thick  lower  layer  was  developed.   This  model  was  used 
to  determine  the  disturbance  structure  in  y  for  linear  inter- 
actions with  a  constant  mean  wind  and  nonlinear  interactions 
with  a  varying  mean  wind.   Simple  heating  to  maintain  the 
mean  wind  and  friction  were  included  in  the  nonlinear  case. 

The  linear  experiment  produced  exponential  growth  rates 
in  the  disturbance  field  that  were  independent  of  the  lower 
layer  thickness.  These  waves  were  baroclinically  unstable. 
The  mean  wind  field  showed  large  modifications  for  a  Rossby 
deformation  length  of  75  KM  at  10^  Jupiter  days. 
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I.   INTRODUCTION 

The  meridional  belts  (dark  bands)  and  zones  (light  bands) 
of  Jupiter  have  been  recorded  since  this  planet  was  first 
observed  through  a  telescope.   The  features  found  in  these 
regions  are  semi-permanent  and  have  a  variety  of  rotation 
rates.   Peek  (1958)  assembled  this  information  to  show  the 
variation  of  zonal  velocity  with  Jovian  latitude.   Evidence 
was  found  by  Chapman  (1969)  for  anticyclonic  shear  to  be 
associated  with  the  zones  and  for  cyclonic  shear  to  be 
associated  with  the  belts,  but  the  correlation  was  imperfect. 

Ingersoll  and  Cuzzi  (1969)  have  investigated  the  baro- 
tropic  instability  of  the  observed  motions  at  mid-latitudes 
using  data  from  Peek  (1958).   Geostrophic  balance  and 
symmetric  temperature  difference  between  the  belts  and  zones 
were  assumed.   The  validity  of  the  assumptions  was  verified 
by  the  excellent  agreement  between  predicted  and  observed 
zonal  velocities.   It  was  also  shown  that  the  relative 
vorticity  gradient   U"  approaches,  but  does  not  exceed  the 
planetary  vorticity  gradient  B.   This  is  the  necessary 
condition  for  barotrophic  instability.   Thus  barotrophic 
instability  probably  limits  either  the  belt  wave  number 
or  the  amplitude  of  the  zonal  flow.   It  was  further  shown 
that  the  zones  are  warmer  than  the  belts. 

Stone  (1967)  applied  the  baroclinic  stability  theory 
to  the  dynamics  of  the  Jovian  atmosphere.   He  tested  the 


hypothesis  that  the  zonal  motions  in  Jupiter's  atmosphere 
are  thermal  winds  and  that  the  latitudinal  cloud  bands  are 
caused  by  baroclinic  instabilities  under  non-geostrophic 
conditions.   Stone  (1966a)  showed  that  at  least  three 
different  kinds  of  baroclinic  instability  can  occur  when 
de/dZ  >  0.   The  kind  of  instability  that  dominates  depends 
upon  the  value  of  the  non-dimensional  Richardson  number 


Ri  =  f!f/(§>2>  (1-15 


where  g  is  the  acceleration  of  gravity,  6  is  the  potential 

1    J\  Q 

temperature,  —  ~y  is  a  representative  vertical  stability, 

3U 
and  7T7T  a  characteristic  vertical  shear  of  the  wind.   The 

o  L 

Richardson  number  is  a  characteristic  ratio  of  work  done 
against  gravitational  stability  to  energy  transferred  from 
mean  to  turbulent  motion.   He  assumed  that  the  observed 
winds  on  Jupiter  are  due  to  symmetrical  instabilities  with 
no  longitudinal  variations.   This  type  of  instability  domin- 
ates in  the  range   0.25  <  Ri  <  0.95-   Allowance  for  a  deep 
atmosphere  was  made  by  not  using  the  Boussinesq  approximation 
for  compressibility.   It  was  shown  that  in  a  deep  atmosphere 
symmetric  instabilities  form  in  the  specified  range  of  Ri 
and  that  these  instabilities  may  extend  over  many  scale 
heights.   Thus  the  zonal  winds  and  the  latitudinal  cloud 
bands  may  be  explained  qualitatively  by  assuming  the  presence 
of  north-sourth  temperature  gradients  caused  by  differential 
solar  heating.   These  gradients  will  drive  thermal  winds  in 


the  zonal  direction,  and  under  appropriate  conditions, 
symmetric  instabilities  can  grow,  generating  latitudinal 
cloud  bands. 

Gierasch,  Ingersoll,  and  Williams  (1973)  demonstrated 
that  a  cloudy  planetary  atmosphere  at  rest  is  unstable  to 
disturbances  of  large  horizontal  scale.   The  energy  source 
for  this  instability  is  the  change  in  radiative  heat  flux 
associated  with  a  vertical  displacement  near  the  emitting 
level.   The  important  feature  of  the  model  was  that  the  net 
heating  rate  Q  is  proportional  to  the  vertical  displacement 
6Z  of  the  cloud  from  its  equilibrium  position.   Quasi-steady 
geostrophic  motions  or  short  period  oscillations,  representing 
inertia-gravity  waves  may  become  unstable  in  the  presence 
of  radiative  heating.   Gierasch  and  Goady  (1969)  found  that 
the  radiative  time  constants  of  Jupiter  are  longer  than  the 
dynamical  time  constants  in  the  earth's  atmosphere;  there- 
fore, the  geostrophic  mode  was  applied,  since  the  instability 
can  be  quasi-geostrophic.   Gierasch,  Ingersoll,  and  Williams 
(1973)  derived  a  value  of  19,000  KM  for  the  y-scale  wavelength 
at  23°  latitude  which  corresponds  to  the  limiting  amplitude 
for  an  axisymmetric  geostrophic  instability  where  the  wave- 
number  in  the  x-direction  k  =  0.   This  y-length  is  essen- 
tially  equal  to  the  observed  belt-zone  pair  at  this  latitude. 
They  also  computed  the  value  of  the  radiative  time  constant 
X  ""  .   For  an  isothermal  atmosphere  X  was  approximately 

equal  tc  ISO   years.   It  was  also  observed  that  even  if  the 
lapse  rate  differed  from  the  adiabatic  lapse  rate  by  only 


1% ,    the  radiative  time  constant  would  be  about  1.6  years. 
These  results  are  consistent  with  the  observed  lifetimes 
of  atmospheric  features  on  Jupiter. 

Gierasch  (1973)  studied  the  radiative  instability  of 
a  cloud  deck  on  an  equatorial  3-plane  as  discussed  by 
Gierasch,  Ingersoll,  and  Williams  (1973).   The  model  was 
based  on  the  assumption  that  a  precipitating  condensable 
atmospheric  constituent  can  be  concentrated  in  the  atmos- 
phere by  vertical  motion.   If  the  condensing  constituent  is 
radiatively  important,  it  can  cause  radiative  warming  at 
the  locations  where  it  is  concentrated.   As  a  result,  the 
heating  is  related  to  vertical  motion  and  reinforcement  is 
possible.   They  considered  two  values  of  the  Rossby 
deformation  length 


a-V^T  • 


Ri  =  Y  &_ii    ,  (1.2) 

where  g'  is  reduced  gravity,  f  is  the  Coriolis  parameter, 
and  H  is  a  characteristic  height.   Static  stability  esti- 
mates by  Lewis  (1969)  and  Stone  (1972)  gave  characteristic 
lengths  of  785  KM  and  78.5  KM,  respectively.   It  was  con- 
cluded that  the  radiative  instability  mechanism  driven  by 
the  dynamical  concentration  of  a  condensable  constituent, 
produces  instabilities  with  the  wavelength  and  zonal  symmetry 
of  Jupiter's  bands.   Nonsymmetric  instabilities  with  a 
shorter  wavelength  also  exist.   Finally,  the  structure  of 


the  modes  is  in  agreement  with  Ingersoll  and  Cuzzi's  (1969) 
conclusion  that  cloudy  zones  are  warmer  than  belts. 

A  three-level  quasi-geostrophic  model  incorporating  a 
deep  lower  layer,  as  suggested  by  Gierasch  (1973),  was 
developed  based  on  the  Charney  and  Phillips  (1953)  multi- 
level model.   The  model  was  modified  by  implementing  Lorenz's 
(I960)  technique  of  replacing  the  dependent  variable  with 
Fourier  coefficients. 

The  objective  of  this  study  was  to  define  the  linear 
disturbance  structure  in  y  and  to  see  if  the  belt-zone 
wind  growth  rate  due  to  quasi-geostrophic  baroclinic 
instability  decreases  with  greater  lower  layer  depth.   Non- 
linear interactions  with  the  belt-zone  flow  pattern  were 
then  studied  incorporating  simple  heating,  to  maintain 
the  zonal  winds,  and  friction.   These  objectives  were 
evaluated  using  two  values  for  the  Rossby  deformation 
length. 


II.   THE  MODEL  AND  FORECAST  EQUATIONS 

A.   DESCRIPTION 

The  atmospheric  model  for  these  experiments  was  based 
on  the  quasi-geostrophic  equations.   The  model  was  constructed 
on  a  3-plane  with  a  central  latitude  of  ^5  degrees.   Since 
the  large-scale  motions  in  the  extra-tropical  latitudes  of 
Jupiter  are  quasi-geostrophic,  the  domain  of  the  model  was 
restricted  to  the  middle  latitudes.   Phillips  (1963)  dis- 
cussed geostrophic  motions  with  approximately  uniform 
potential  vorticity.   These  motions  have  characteristic 
wavelengths  of  about  4,000  KM  on  earth,  or  the  order  of 
the  Rossby  deformation  length.   He  used  an  expansion  in 
the  Rossby  number,  R  =  U/fL,  where  U  is  a  characteristic 
velocity,  f  is  the  Coriolis  parameter,  and  L  a  character- 
istic length.   This  expansion  showed  that  the  $-plane 
approximation  is  justified  for  geostrophic  motions  of  this 
type.   Gierasch,  Ingersoll,  and  Williams  (1973)  have  shown 
that  large-scale  bands  have  a  characteristic  scale  of  15,000 
KM  in  the  atmosphere  of  Jupiter.   It  was  further  shown  that 

the  Rossby  number  R  is  much  less  than  one.   If  the  rota- 
J  o 

tion  is  sufficiently  fast  such  that  R  is  of  order  one  or 

o 

less,  Coriolis  forces  become  important  and  must  be  considered 
•in  the  equations  of  motion.   The  additional  conditions  of 
RT/a,  where  a  is  the  planetary  radius,  being  small,  and  the 
large-scale  waves  being  greater  than  or  equal  to  RL  are 
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satisfied;  therefore,  the  geostrophic  assumption  is  valid. 
The  geophysical  constants  for  earth  and  Jupiter  are  summa- 
rized in  the  appendix. 

A  three-level  model  was  developed  so  that  heating  and 
friction  could  be  included  and  still  satisfy  the  lower 
boundary  condition.   Since  the  interior  of  the  planet  is 
not  known,  a  deep  lower  layer  is  used  to  remove  the  two 
upper  layers  from  the  lower  boundary.   The  solution  function 
and  its  derivatives  were  forced  to  vanish  at  the  boundaries. 
The  Jovian  atmosphere  was  divided  into  six  layers  numbered 
1  through  6  as  shown  in  Figure  1.   The  upper  layers,  0-4, 
have  thickness  AP/2 ;  and  the  lower  layers,  5  and  6,  have 
thickness  P/2,  where  P>>AP  in  this  application.   The  thick- 
ness between  levels  3  and  5  is  defined  to  be  n  =  (AP+P)/2. 
In  general  these  thickness  relationships  are  arbitrary. 

B.   METHOD  OF  SOLUTION 

The  quasi-geostrophic  vorticity  equation  on  a  3-plane  is 

f£+  v.vcc+e0y).-  f0  §£-  bv2c  ,  (2.D 


where      c   =   ■£-  V2<})    ,        V  =   —-  kxV<f>    ,  (2.2) 

o  o 


and  the   thermodynamic   equation   is 


3      3<fc    ,    ,,.  „  36    ,      _        1    /3<t)*        3<K  ro    ?\ 

It    DP   +   V'V    9P    +    Wa    =    T   (aP"  "   W)       ■  (2*3) 
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FIGURE    1    -  The   Three-Level   Model 
Variables  are  evaluated   at   the   layers    shown. 
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dp 
where  co  =  tt  and  a  is  defined  as  the  static  stability.   The 
at 

quantity  t  is  the  radiative  time  constant  and  — ^—  is  an 

d? 

equilibrium  temperature  which  is  dependent  upon  radiation, 
convection,  and  other  small  scale  parameters.   Equation  (2.1) 
was  applied  at  levels  1,  3,  and  5  which  yielded 

^  +  VV(c1+B0y>  ~  f0  pppf  =  BV \     >  (2^} 

■^  +  V0.V(5,+3J)    -    f      p^-p^  -    BV^Co       ,  (2.5) 


3t  v3      V4,3   Mc"  o  P^-P2       ~¥    s3 

^  +  VV( VBoy)  -  fo  pp^  =  bv%    •  (2-6) 

Equation  (2.3)  applied  to  levels  2  and  4  yielded 

(p-Zp-)  +  V2-V  (p-zp-)  +  a2o)2  -  ~   [(Fp-^  -  (^)2]   , 


at  ^3^        .3  u 


(2.7) 


3t  (P^pf}  +  VV  (P7^pt)  +  G^4  "  7  C(8P  ^  "  (3pV   • 

(2.8) 


Equations  (2.7)  and  (2.8)  were  solved  for  w2  and  co^, 
respectively,  and  substituted  with  Equation  (2.2)  into 
Equations  (2. A),  (2.5),  and  (2.6).   These  operations  yielded 
the  potential  vorticity  equations 
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~   [V201  +  k2(c})3-(})1)]    =    -V1-V(V2(j)1  +  fo3oy)    -    K2V2.V(cf)3-cf)1) 

+    BV%1    +   -i[34,2-((j)3-(j)i)]    =    ci    ,     (2.9) 


|^  [V2(^3+K1|((|)5-cf)3)-K2((})3-({)1)]  =   -V3.v(v24>3+fo30y) 


-    k^V^- V(<f>5-4>3)    +    k2V2-V(<J)3-<J)1) 


+  bv%3  +  JL[a<(,J-((|,5-<(>3)] 


^[34,*- C^-^)]  =  c. 


(2.10) 


|^  [v2(j)5-K6((t)5-4»3)]  =  -v5-v(v24»5+f030y)  +  k6v4-v(*5-*3) 


+    BV%5    -   —  [94>J-U5-$3)]    =    C3  ,  (2.11) 


where 


2 


'2         At)2  5 

AP   a2 


(2.12) 


and 


Kh    = 


APttcj 


=    6k, 


Pa6" 


=     £K, 


(2.13) 


(2.14) 
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A  set  of  three  linearly  independent  eigenvector  equa- 
tions were  formed  to  solve  the  equation  set  (2.9),  (2.10), 
and  (2.11)  after  a  method  by  Charney  and  Phillips  (1953). 
If 

(2.15) 


then  the  general  form  of  the  potential  vorticity  equation 

is 

v  at  +  h  K2  at  =  5  >  (2-l6) 


where 


and  A  is  the  appropriate  matrix.   If 


(2.17) 


is  an  eigenvector  of  A  and  X  is  its  corresponding  eigenvalue, 
then 

Ay  =   Xy  (2.18) 
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and 


-1 

1 

0 

1 

-(1+6) 

£ 

0 

6 

-E 

A      =  1  -(1+6)  e  ."  (2.19) 


Now   if 

•      =    HliLl    +   H2ii2    +   H3y_3         ,  (2.20) 

and  C      =   d.y,    +   d^u^   +   cUjj~         ,  (2.21) 

are  substituted  into  Equation  (2.16)  and  Equation  (2.18)  is 
used  to  determine  the  eigenvalues,  then  the  coefficients  of 
each  eigenvector  can  be  equated  to  give 


ft  C(V2+AiK2)Hi]   =  d±         ,  (2.22) 


where   i  =  1,2,3. 

The  following  eigenvalues  as  specified  by  Equation  (2.18) 
were  obtained, 


*1   =   0 


X2   =   "  I  ^2  +  6  +  e)  +  |  [(e  +  6)2  +  ^(l-e)]i5     \     .      (2.23) 


X3   =  '-  |  (2  +  6+e)  -  \   C(e+6)2+Ml-0]^ 
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The  eigenvector  equation  set  was  formed  from  the  poten- 
tial vorticity  equations  (2.9),  (2.10),  and  (2.11)  and 
Equation  (2.22).   These  equations  are 


|jT  [V2H1]  =  -V-[6oykxVH1+V1V2(f)1+V3V2(|)3+V5V2|<|)5] 

+  BV2H1  +  Q     ,  •  (2.24) 


|t    W   +X2k2)H2]    =    -V-{BoykxVH2+V1V    4>1+32V  V24)3+y2V5V2c^5 

-    X2K2[V2    ^3_(l)i)-Y2V^(*5~*3^}    +    ByilH3   +    Q      ' 

(2.25) 


~  [(V2+X3k2)H31  =  -V(BoykxVH3+V1V2(f)1+^3V3V2(j)3+Y3V5V2({)5 

-  ^3k2Cv2  (^3-4>1)-Y3Vi4  (4)5-<J>3)l>  +  BvV5  +  Q   , 

(2.26) 

where  Q  *  -~  [  (1-3 )  (3<f»2-3<f>2 )  +  OB-ye)  O^J-B^)]  .    (2.27) 


An  assumed  solution  of  the  form 

H  =  H(y,t)  +  H  cos  kx  +  H   sin  kx   ,  (2.28) 

c  s 

where  H  is  the  zonal  mean  solution  amplitude,  H  is  the 
Fourier  cosine  eddy  amplitude,  and  H  is  the  Fourier  sine 
eddy  amplitude,  was  introduced  into  Equations  (2.24),  (2.25), 
and  (2.26)  as  suggested  by  Lorenz  (I960).   The  equations 
were  quasi-linearized  by  eliminating  all  terms  of  higher 
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wave  number  than  k.   This  technique  permits  full  interaction 
of  the  disturbance  in  x  with  the  mean  flow,  but  neglects 
interaction  with  the  other  waves  in  x. 

Equations  (2.24),  (2.25),  and  (2.26)  were  separated 
into  components  independent  of  x,  coefficients  of  cos  kx, 
and  coefficients  of  sin  kx.   The  resulting  equation  set 
(2.29)  through  (2.^3)  constitute  the  prediction  equations 
for  this  model.   The  primed  quantities  represent  the  partial 
derivative  with  respect  to  y. 

The  terms  independent  of  x  in  (2.24)  give: 


£  t!  V  til  111  H  t  Ml 

9t   H      =    -   2f^  [(*lc*ls-({)lshc)-((t,lshc-<{,lc(t)ls) 


tii  in  ii       i         n       i 

+  (*3cc()3s~<i)3s(f)3c)-U3scf)3c-*3c*3s) 


r  ill  IM  £  ii  I  ii         it 

+  e   ^c^s-^s^c^  z    (*5b*5c"*5c*5b)] 


+   F  +   Q         ,  (2.29) 


where  F  =   BH  (2.30) 


and  Q  =   -j-  [(l-3)(34>*-^2)  +  (6B~Ye)(ac|)*~3^l):]    *         (2-3D 
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The  cos  kx  terms  in  (2.24)  give: 


ft    (^2   "   k2)Hc    '    "   r^oVo   +    ^l^hhs-h^ 

+    ^3M+k2*>3s-*3*3s) 


+    |   (^M+k2?^5s-^^s)]    +   Pc    +    Qc    , 

(2.32) 
where  F      =   B(H""-2k2H"+kV    )         ,  (2.33) 


and      Qc  =  —  [(l-B)04'2-3*2c)  +  ((53-Yc)(34»*-8({)4c)]  .(2.34) 
The  sin  kx  terms  in  (2.24)  give: 

|_  (iL  .  k2)Hs  =  -  f  [BQf0Hc  +  (*;%10+k2*Xc-*Xc) 
3y  o 

+  (^%3c+k2^^c-?^3c) 


•    + 1  (*;%5c+k2?;^5c-^^c)]  +  fs  +  qs, 

(2.35) 

where    F   =  B(H""-2k2H"+k  H  )    ,  (2.36) 

s       s       s     s     ' 

and      Qg  =  -£   C(l-3)04>2-3*2s)  +  (63-Ye)0<|>IJ-3<I)lls)]  .(2.37) 
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The  terms  independent  of  x  in  (2.25)  give: 


ft  (^2  +  V2)H  =  "  if  t<*i>i8-*i>io)-(*lB*io"*lC*i8) 


+Y2[(*5c*5s-*5s*5c)-(,|'5s*5c-','5o*5s):l 


X2K2C(,|,3o*ls-*3s,t'lc)-(*lo'|,3S-,|'lS*3c) 


-Y2(*5c*3B-*5s*3o)+Y2(*3c*5s-*3s*5c):i1 


+  F  +  Q     ,  (2.38) 


The   cos   kx  terras   in    (2.25)    give 


2 

If   <^7  -k2  +  X0K?)H  =    -   -£-  {8fH      +    ($'%      +k2^Y    -?nV  _) 

dt        ~     2                     2     2C                  f  O    O     S  1        IS              1     IS        1     IS 

3y  o 


+Y2(^%5s+k2^<f)5s-^^s) 


+    X2K2[(^^3s^^ls)+y2(^*3s^>5s)]}    +   Fc    +    Qc, 


(2.39) 
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The    sin   kx  terms    in    (2.25)    give: 


2 

It    (^2    -*2+V2)Hs    "    "   £   <VoHc    +    (*i"*lc+k^>ic^>ic) 


+Y2(^5%5c+k2c()^5c-^^c) 


+  X2K2[(*>3c^3*lc)    +  ^2(*^3c^>5c)]}   +   Fs   +   Qs  > 

(2.40) 

The  terms  independent  of  x  in  (2.26)  give: 


2 

|_  (»   +x  K  )ff  =  -  jk-  {(♦i>l8-Olc)-(*l8*ic-*le*i.) 

9y  o 


^3C(*5c*5S-*5>5c)-(*5s*5c-+5c*5S)] 


-X3K2[(*3o*ls-*3s*lo)-(*lC*3s-*ls*3c) 


-V*5c*3s-*5s*3c)+Y3(*3c*5s-*3s*5c)]> 


+  F  +  Q   ,  (2.11) 
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The  cos  kx  terms  in  (2.26)  give: 


2 

If  (Z2  -k2+x3K2)Ho  =  -  f  te0foHs  +  (*r*is+k2*i*is-*i*is) 

ay  O 


+X3ic2C(<|)1({)3s-(()3(J)1s)+y3((|>5<J)3s-<{>3<()5s)]}  +  Fc  +  Q 


c  » 
2.42) 


The  sin  kx  terms  in  (2.26)  give: 


2 

_  (i_  -k2+x3.2)Hs  -  -  Jl  {6ofoHc  +  (♦i,*lc+k2y^lc.y^;c) 

+  e3(^%3c+k2^cf)3c-*3cD3c) 


+Y3(*5,*5c+k2*5*5c"*5*5c) 


+  X3K2[(^^3c^3$lc)+y3(*^3C^3^c)]}  +  Fs  +  Qs  . 

(2.43) 

A  method  developed  by  Richtmeyer  (1957)  was  used  to 
solve  for  the  tendency,  Tjr-  ,  at  each  time  step.   The  equa- 
tions were  integrated  in  time  by  a  marching  process.   All 
finite  differences  were  centered  in  time  or  space.   A 
forward  time  step  was  used  to  start  the  marching  process. 
A  finite  difference  scheme  developed  by  Matsuno  (1966)  was 


22 


used  every  50  time  steps.   This  is  a  two  step  iteration  to 
simulate  the  backward  difference  method.   This  scheme  was 
applied  in  this  situation  to  suppress  the  computational 
mode  that  would  arise  from  using  the  centered  time  differ- 
encing exclusively.   A  time  step  of  30  minutes  was  used. 
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III.   RESULTS 

Two  experiments  were  conducted  under  varied  conditions. 
The  first  experiment  was  designed  to  yield  the  growth  rate 
of  the  most  unstable  wave  in  the  disturbance  fields  for 
Rossby  deofrmation  lengths  of  75  KM  and  750  KM.   This  range 
of  values  for  the  Rossby  deformation  length  was  suggested 
by  Gierasch  (1973).   The  upper  two  levels  each  had  a  pressure 
differential  of  50  CB,  while  thickness  of  the  lowest  level 
assumed  values  of  50,  500,  and  5000  CB,  respectively.   The 
y-wavelength  of  the  mean  flow  was  15,000  KM  for  all  cases. 
The  x-wavelength  varied  from  1,000  KM  to  15,000  KM  for 
R   =  750  KM  and  from  300  KM  to  3000  KM  for  RL  =  75  KM. 
k?  as  defined  by  Equation  (2.12)  was  scaled  to  be 

<P  =   -ip  .  (3-D 

RL 

The  static  stability  at  level  2,  a-,  was  obtained  from 
Equation  (3.1)  and  the  static  stability  at  level  4,  a^, 
was  found  by  the  relation 

P22 

ar  a2  72  '  (3-2) 

where  P?  and  P^  are  the  absolute  pressures  at  levels  2  and 
4,  respectively.   This  ratio  was  a  constant  equal  to  0.25 
for  all  experiments.   The  mean  wind  U  was  initialized  as 
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a  sinusoidal  function  of  y  with  an  amplitude  of  20.5 

meters  per  second  in  level  1.   This  value  was  obtained 

2 
from  U   =  3/k   which  is  the  limiting  value  for  barotrophic 

stability  as  suggested  by  Ingersoll  and  Cuzzi  (1969).   The 

mean  wind  was  zero  in  levels  3  and  5.   The  mean  field  was 

held  constant,  which  makes  the  equations  linear,  throughout 

the  first  experiment.   The  disturbance  field  cos  kx  contained 

an  initial  sinusoidal  variation  with  the  same  y-scale  as 

the  mean  flow  and  an  amplitude  an  order  of  magnitude  smaller 

than  the  mean  field.   The  initial  sin  kx  field  was  entirely 

described  by  the  cos  kx  terms.   There  were  no  heating  or 

friction  terms  included. 

The  forecasts  were  carried  forward  until  the  growth 

rates  became  exponential.   Growth  rates  of  the  disturbance 

were  calculated  by 


1  3A  _  M  _    1      ,   A2 

A  It  ~  M  "  t2  -  tx  ln  A^  ,  (3.3) 


where  M  is  the  exponential  growth  rate,  A  is  the  amplitude 
of  the  disturbance,  and  t«  and  t,  are  the  times  at  which 
the  disturbance  amplitudes  were  A~  and  A,,  respectively. 
The  maximum  growth  rate  for  R,  =  750  KM  occurred  at  an 
x-wavelength  of  6,000  KM  for  all  lower  layer  thicknesses 
(Figure  2).   For  lower  layer  thicknesses  of  50  and  500  CB 
and  RT  =  75  KM  the  maximum  growth  rate  was  at  an  x-wavelength 
of  600  KM.   The  5000  CB  case  showed  maximum  growth  rate  at 
an  x-wavelength  of  500  KM  (Figure  3). 
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FIGURE   2   -   Growth  Rate   of   Y-Disturbance  Wave   for 
RL=750KM 
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The  phase  and  amplitude  relationships  of  the  disturbance 
field  were  determined  by 


<J>      cos   kx   +    <J>      sin   kx   =    C    cos(kx   -   6), 

(3.*) 

where  <j>  and  A  are  the  geopotentials  in  the  cos  kx  and 
sin  kx  fields,  respectively.   Equation  (3.4)  was  expanded 
to  yield  the  following  expressions  for  the  amplitude  C 
and  the  phase  angle  <5 , 


6  =  arctan  (^)  ,  (3.5) 
c 

C  =  j —  =  — : t —  .                   (3.6) 

cos  6    sin  6  VJ  ' 


The  amplitude  and  phase  relationships  for  A..  ,  A, .  >    and 
<j)j-  were  computed  at  the  maximum  growth  rate  for  each  case. 
The  results  are  shown  in  Figures  4  through  15.   All  ampli- 
tude profiles  showed  greater  relative  growth  in  the  northern 
wave;  however,  the  effect  was  more  pronounced  at  RT  =  750 

Li 

KM.   This  was  due  to  the  presence  of  the  planetary  vorticity 
gradient  3.   The  relative  amplitudes  of  A   A   and  A 
were  the  same  order  of  magnitude  for  the  lower  level  thick- 
ness P  =  50  CB  (Figures  4  and  7).   The  amplitude  of  A 
was  an  order  of  magnitude  smaller  than  A  or  A  for 
P  ^   500  CB  (Figures  5  and  8).   At  P  =  5000  CB  (Figures  6 
and  9)  the  amplitude  of  Ar  was  negligible  compared  to  A 
or  A  .   These  characteristics  were  observed  at  both  values 
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FIGURE    k    - 


Relative  Amplitude  of  <{>,,  <K,  and  <}>,- 
as  a  Function  of  Y  for  AP=50  CB, 
R  =  75  KM,  and  x-wavelength  of  600  KM, 
Y-Axis  has  been  shifted  to  the  north. 
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FIGURE  5  -  Relative  Amplitude  of  <J>. 
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FIGURE  6  -  Relative  Amplitude  of  4>1  , 


(J)^,  and  <}>,-. 


as  a  Function  of  Y  for  AP  =  50  CB, 
P  =  5000  CB,  RL  =  75  KM,  and 

x-wavelength  of  600  KM. 
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FIGURE  7  -  Relative  Amplitude  of  <J>. 
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as  a  Function  of  Y  for  AP  =  50  CB, 
P  =  50  CB,  RL  =  750  KM,  and 
x-wavelength  of  6000  KM. 
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-   Relative   Amplitude    of   <K  ,    cf)^,    and   <{>,- 

as   a  Function   of  Y   for   AP   =   50   CB, 
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FIGURE  10  • 

-  Phase  Angles  of  <f>-.,  <f>_,  and  cfv 
Function  of  Y  for  AP  =  50  CB, 
P  =  50  CB,  RL  =  750  KM,  and 
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FIGURE    11    ■ 

-   Phase   Angles   of   4>n  ,    <J>_,    and   §c 
Function   of  Y   for   AP   =    50   CB, 
P   =    500    CB,    RL   =    750   KM,    and 
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FIGURE  13  - 


Phase  Angles  of  <J>,,  <$>*,    and  <{>,-  as  a 
Function  of  Y  for  AP  =  50  CB, 
P  =  50  CB,  RL  =  75  KM,  and 
x-wavelength  of  600  KM. 
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FIGURE  14  - 


Phase  Angles  of  4>,,  4>_,  and  <{>,-  as  a 
Function  of  Y  for  AP  =  50  CB, 
P  =  500  CB,  RL  =  75  KM,  and 
x-wavelength  of  600  KM. 
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FIGURE  15  - 


Phase  Angles  of  <K,  <J>_,  and  $>,-   as  a 
Function  of  Y  for  AP  =  50  CB, 
P  =  5000  CB,  RL  =  75  KM,  and 
x-wavelength  of  600  KM. 
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of  the  Rossby  deformation  length.   The  phase  angles  for 
each  level  showed  that  all  of  the  waves  were  barotropically 
stable,  since  the  phase  lag  in  y  was  in  the  same  sense  as 
the  horizontal  wind  shear.   Also,  the  vertical  tilt  was  in 
the  opposite  sense  to  the  vertical  shear  as  is  to  be  ex- 
pected for  baroclinic  instability.   Note  that  the  y-scale 
of  the  disturbance  is  much  smaller  for  PL.  =  75  KM  compared 
to  RL  =  750  KM.   This  result  was  predicted  by  Stone  (1969) 
who  showed  that  in  a  slowly  varying  wind  field  the  y-scale 
of  baroclinic  waves  would  be  of  the  order  of  the  Rossby 
deformation  length.   It  was  observed  that  the  amplitude 
in  the  disturbance  fields  for  RT  =  75  KM  was  contained  in 
two  narrow  regions  centered  about  2600  KM  and  7500  KM. 
Friction  and  heating  were  added  to  the  model  in  the 

second  experiment  to  study  the  nonlinear  effects.   The 

%  fh  % 
equilibrium  temperature  profile  —-   was  obtained  from 

the  initial  value  of  -^   .   This  has  the  effect  of  main- 
taining the  mean  wind,  since  the  heating  will  tend  to 
restore  the  mean  temperature  field  to  its  initial  conditions 
when  the  disturbance  amplitude  is  sufficient  to  cause  a 
departure.   Heating  was  present  in  the  mean  field  only. 
The  radiative  time  constant  x  was  10  Jupiter  days.   The 
friction  term  was  present  in  all  fields.   The  friction 
coefficient  B  was  equal  to  1.0  x  10   and  5.0  x  10^  for 
Rossby  deformation  lengths  of  750  KM  and  75  KM,  respectively 
Since  the  grid  mesh  distance  was  150  KM  and  the  x-wavelength 
was  600  KM,  the  grid  may  have  difficulty  resolving  the  wave; 
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therefore,  a  larger  value  of  the  friction  coefficient  was 

needed  to  maintain  stability  in  the  model.   The  mean  field 

was  allowed  full  interaction  with  the  disturbance  fields. 

Two  cases  corresponding  to  the  Rossby  deformation  lengths 

were  performed  with  a  lower  layer  thickness  of  500  CB. 

All  other  initial  conditions  were  as  described  in  the  linear 

experiment. 

The  model  was  integrated  for  104  Jupiter  days  with  a 
Rossby  deformation  length  of  75  KM.   At  this  time  the  model 
became  unstable.   The  corresponding  mean  zonal  wind  profile 
and  disturbance  amplitudes  and  phase  relationships  are 
shown  in  Figures  16  through  18.   The  maximum  wind  speeds 
were  25.9  meters  per  second  for  the  northern  easterlies  and 
26.3  meters  per  second  for  the  southern  westerlies.   Zonal 
velocities  also  developed  in  layers  3  and  5.   The  relative 
minima  associated  with  the  maximum  wind  is  the  result  of 
the  amplitude  of  the  disturbance  fields.   These  amplitudes 
are  similar  to  the  linear  case  except  the  peaks  are  now 
bimodal  and  the  profile  has  a  broader  base,  indicating  the 
disturbance  energy  has  been  distributed  over  a  larger 
region.   The  phase  angles  showed  baroclinic  instability 
in  the  vicinity  of  the  maximum  mean  zonal  wind  in  the 
north  and  the  south,  since  the  sense  of  the  vertical  wind 
shear  was  opposite  the  tilt  of  the  phase. 

Because  the  750  KM  Rossby  deformation  length  case  has 
a  much  slower  disturbance  growth  rate,  the  experiment  was 
run  for  536  Jupiter  days.   The  mean  zonal  wind  profile  is 
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FIGURE   16   -  Mean   Zonal   Wind    for   AP   =    50   CB, 
P   =    500   CB,    RL   =    75   KM,    and 
x-wavelength   of   600   KM   for 
Jupiter  Day   104. 
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FIGURE  17  -  Relative  Amplitude  of  <J>. 


4>o>  and  <f)[- 


as  a  Function  of  Y  for  AP  =  50  CB, 
P  =  500  CB,  RL  =  75  KM,  and 
x-wavelength  of  600  KM  for 
Jupiter  Day  10*1. 
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FIGURE    18    - 


Phase   Angles   of   <J>-  ,    <}>_,    and   <$>c   as   a 

J-    D  D 

Function  of  Y  for  AP  =  50  CB, 
P  =  500  CB,  RL  =  75  KM,  and 
x-wavelength  of  600  KM  for 
Jupiter  Day  104. 
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shown  in  Figure  22.   Relatively  small  changes  have  occurred 
compared  with  the  other  nonlinear  case.   The  mean  wind 
speed  is  20.7  meters  per  second  in  level  1.   Small  amplitude 
winds  have  developed  in  levels  3  and  5;  however,  no  appre- 
ciable wind  is  observed  in  level  5  in  the  south.   Figure  21 
shows  the  amplitude  of  the  most  unstable  wave  in  the  distur- 
bance field  as  a  function  of  time.   The  growth  rate  is 
exponential  with  value 


1   3A        a  on      n  n~6       -1 

■r   rr  =   0. 21  x  10    sec 


thus  the  field  has  not  reached  quasi-equilibrium.   The 
disturbance  amplitudes  shown  in  Figure  19  indicate  the 
effects  of  the  planetary  vorticity  gradient  3  have  been 
reduced.   The  northern  disturbance  has  a  relatively  greater 
amplitude  in  all  three  levels.   The  relative  amplitude  of 
<J>   is  an  order  of  magnitude  smaller  as  observed  in  the 
linear  case. 

One  inconsistency  was  present  in  this  experiment  which 
has  not  been  resolved.   The  growth  rates  observed  in  the 
linear  mode  were  an  order  of  magnitude  larger  than  the 
nonlinear  mode.   This  was  not  a  physical  result,  but  an 
error  in  the  program. 
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FIGURE  19 

-  Relative  Amplitude  of  <{>,  ,  $_,  and 
as  a  Function  of  Y  for  AP  =  50  CB3 
P  =  500  CB,  RL  =  750  KM,  and 
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FIGUR 

E  20  -  Phase  Angles  of  4>,,  <K,  and  cfv 
Function  of  Y  for  AP  =  50  CB, 
P  =  500  CB,  RL  =  750  KM,  and 
x-wavelength  of  6000  KM  for 
Jupiter  Day  536. 
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IV.   CONCLUSIONS 

A  three-level  quasi-geostrophic  atmospheric  model  has 
been  examined  on  a  3-plane  with  a.  single  wave  in  the  x- 
direction  and  a  deep  lower  layer  of  varying  thickness.   The 
disturbance  structure  in  y  was  investigated  to  determine 
the  linear  interactions  with  a  constant  mean  wind  and  for 
nonlinear  interactions  with  a  varying  mean  wind.   Simple 
heating  to  maintain  the  mean  wind  and  friction  were  included 
in  the  nonlinear  case.   Two  values  of  the  Rossby  deformation 
length  were  included  in  the  experiments. 

The  linear  experiment  showed  that  the  growth  rates  in 
the  disturbance  field  are  largely  independent  of  the  lower 
layer  depth  (Figure  2).   These  growth  rates  are  a  function 
of  the  Rossby  deformation  length  and  the  x-wavelength.   The 
waves  were  barotropically  damped  in  this  case,  since  the 
relative  vorticity  gradient  was  less  than  or  equal  to  the 
planetary  vorticity  gradient.   This  was  obvious  since  the 
phase  tilt  was  in  the  same  sense  as  the  horizontal  wind 
shear. 

The  heating  function  was  not  important  in  the  nonlinear 
experiment,  since  the  model  never  reached  a  long  term  steady 
state.   The  disturbance  fields  were  barotropically  damped. 
In  the  cases  studied  the  zonal  disturbances  are  baroclinically 
unstable.   If  the  actual  wind  structure  on  Jupiter  is  baro- 
clinically stable  as  suggested  by  Ingersoll  and  Cuzzi  (1969), 
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then  this  could  be  accomplished  by  distributing  the  mean 
change  through  great  depth.   This  would  give  a  very  small 
shear  which  could  be  stabilized  by  3.   The  growth  rates  were 
expected  to  decrease  with  increasing  lower  layer  depths, 
but  in  fact  they  did  not.   Thus  it  may  be  that  the  lower 
layer  was  not  properly  formulated  by  this  model. 

Further  investigation  of  the  disturbance  growth  rates 
is  required  to  resolve  the  discrepancy  between  the  linear 
and  nonlinear  cases.   It  is  also  recommended  that  the  grid 
mesh  be  reduced  to  more  realistically  model  the  600  KM 
x-wavelength.   This  should  allow  use  of  a  smaller  friction 
coefficient . 
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APPENDIX 
GEOPHYSICAL  AND  ASTRONOMICAL  CONSTANTS  OF  JUPITER  AND  EARTH 


PARAMETER 


Perihelion 
Distance,  A.U.     8 
(1  A.U.  =  1-.496  x  10  KM) 

Amphelion 
Distance,  A.U. 

Siderial 
Period,  (day) 

Inclination  of 
orbit  to  ecliptic 

Mass,  (GM) 

Density,  (GM/CM3) 

Mean  Radius,  (KM) 

Equatorial 

Gravitational        P 
Acceleration,  (CM/SEC^) 

Albedo 

Orbital 
Eccentricity 

Mean  Temperature 

Solar 

Constant,  (Ly,MIN) 

ft,  (SEC"1) 

Observed  Mean 
Wind,  (M/SEC) 

Rossby  Number,  R 

Rossby 

deformation  Length,  (KM) 


JUPITER 

EARTH 

4.951 

0.9833 

5.455 

1.0167 

10,759.: 

20 

365.256 

3°04' 

23°27' 

1.899  x 

1C 

,30 

5.976  x  1027 

1.33 

5.517 

69,828 

6,371 

2587 

979 

0.73 

0.36 

0.0483 

0.01627 

130°  K 

330°  K 

-0.07 

~2 

17.6  X 

10" 

•5 

7.3  x  10~5 

50 

15 

0.0041 

0.032 

785,  78 

.5 

1,000 
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PARAMETER  JUPITER  EARTH 

Coriolis  Parameter,  (45°)  2.5  x  10   SEC-1  1.0  x  10"   SEC-1 

Planetary  Vorticity  ,  _n  3.6  x  10"12  1.6  x  10-11 
Gradient, (45°), (SEC  1M   ) 

RL/a  0.05  0.63 
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